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The quantum Coulomb glass model describes disordered interacting electrons on the insulating side 
of a metal-insulator transition. By taking quantum fluctuations into account it can describe not only 
the localized limit but also the weakly localized regime. We discuss several possibilities to generalize 
the concept of Anderson localization to interacting electron systems such as the quantum Coulomb 
glass and define criteria for localization. The corresponding physical quantities are calculated by 
numerically exact diagonalization. The results indicate that single-particle excitations close to the 
Fermi energy become more strongly localized under the influence of interaction. 
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I. INTRODUCTION 

Disordered interacting electronic systems such as semi- 
conductor impurity bands or granular metals have at- 
tracted a large amount of theoretical and experimental 
interest fij. Even the presence of disorder or interaction 
alone leads to interesting phenomena, and the compe- 
tition between the two gives rise to a very rich behav- 
ior. The single-electron states in non-interacting systems 
become localized in space for strong enough disorder. 
This phenomenon is called Anderson localization it 
is responsible for the metal-insulator transition in non- 
interacting systems. The generalization of this concept 
to interacting systems is not straightforward since many- 
particle physical properties cannot in general be derived 
from single-particle properties. Consequently, a defini- 
tion of localization in many-particle systems will not be 
unique. In this paper we therefore consider several possi- 
bilities to generalize the concept of Anderson localization 
to interacting systems and discuss their relative merits. 
As an example we then apply the resulting localization 
criteria to the quantum Coulomb glass model of disor- 
dered insulators and calculate the corresponding physical 
quantities by numerically exact diagonalization. 

The paper is organized as follows: In Sec. II we moti- 
vate and define the quantum Coulomb glass model. Sec. 
Ill is devoted to a discussion of Anderson localization in 
interacting systems and the presentation of our results. 
Sec. IV contains some discussions and conclusions. 



II. FROM CLASSICAL TO QUANTUM 
COULOMB GLASS 

Almost the entire current understanding of disordered 
insulators has been obtained from studying the insulat- 
ing limit where the electrons are completely localized 
and can thus be described as classical point charges (we 



will denote this regime as the classical insulating regime 
from now on). Although it is not a quantum mechani- 
cal system, the classical disordered insulator is a com- 
plicated many-body problem. Pollak predicted |s| an 
interaction-induced reduction of the single-particle DOS 
at the Fermi energy in disordered insulators. Later Efros 
and Shklovskii defined Q the generic model of the classi- 
cal insulating regime, the classical Coulomb glass model 
which consists of point charges in a random potential 
which interact via Coulomb interactions. The model is 
defined on a regular hypercubic lattice with N = L d (d is 
the spatial dimensionality) sites occupied by KN (spin- 
less) electrons (0 < K < 1). To ensure charge neutrality 
each lattice site carries a compensating positive charge 
of Ke. The Hamiltonian of the classical Coulomb glass 
reads 



H c i = ^{f i - y)rn + \ y^(nj - K)(nj - K)U K 



(1) 



Here rij is the occupation number of site i and /i is 
the chemical potential. The Coulomb interaction Uij — 
e 2 /rij remains long-ranged since screening breaks down 
in the insulating phase. We set the interaction strength 
of nearest neighbor sites to 1 which fixes the energy 
scale. The random potential values ipi are chosen inde- 
pendently from a box distribution of width 2 Wo and zero 
mean. The physics of the classical Coulomb glass model 
has been investigated by several analytical and numer- 
ical methods and its properties are comparatively well 
understood by now || although the nature of the trans- 
port mechanism is still controversially discussed . One 
of the remarkable features is the power-law gap in the 
zero-temperature single-particle density of states (DOS) 
which is called the Coulomb gap 

Since experiments deep in the insulating regime are 
difficult to carry out most results on disordered insula- 
tors have been obtained from samples not too far away 



1 



from the metal-insulator transition S. Here the (single- 
particle) localization length is still much larger than the 
typical distance between two sites and the description 
of the electrons in terms of classical point charges be- 
comes questionable. In order to investigate the influence 
of finite overlap between the states on the properties 
of the insulating phase we have defined || the quan- 
tum Coulomb glass model, the minimal model of the 
"quantum insulating regime" which accounts for disor- 
der, long-ranged interactions and the quantum nature of 
the electrons. It is obtained from the classical Coulomb 
glass by adding hopping matrix elements of strength t be- 
tween nearest neighbors. The Hamiltonian of the quan- 
tum Coulomb glass reads 

H = -tY,(4 c 3 + c h) + H ch (2) 

<«) 

where cf and Cj are the electron creation and annihilation 
operators at site i, respectively, and the sum runs over 
all pairs of nearest neighbor sites. In the limit t — > the 
model (||) reduces to the classical Coulomb glass, for van- 
ishing Coulomb interaction but finite overlap it reduces 
to the usual Anderson model of localization. 

III. LOCALIZATION IN AN INTERACTING 
SYSTEM 

It has been known for a long time jlO| that disorder 
leads to spatial localization of single-particle states in 
a non-interacting system. Since a system with local- 
ized states at the Fermi level is insulating, the transition 
between delocalized and localized states corresponds to 
a metal-insulator transition. There is, however, much 
experimental and theoretical evidence that a descrip- 
tion of disordered electronic systems in terms of non- 
interacting particles is inadequate. For this reason one 
would like to apply the concept of localization to inter- 
acting systems. Unfortunately, a direct generalization 
of the non-interacting case is not possible since single- 
particle states are not defined in an interacting system 
while the many-particle states always correspond to ex- 
tended charge distributions. The deeper reason for these 
difficulties is, of course, that in a non-interacting sys- 
tem the many-particle properties are completely deter- 
mined by the single-particle properties whereas the same 
is not true for interacting systems. Therefore, in a many- 
particle system one can consider several types of "local- 
ization" (for single-particle or different many-particle ex- 
citations) which are all unrelated a priori. In this section 
we discuss some of these ideas and we also present re- 
sults for the corresponding physical quantities of a two- 
dimensional quantum Coulomb glass. 

From an experimental point of view the most nat- 
ural quantity to consider is probably the conductance 



since it is easily measurable and its behavior determines 
whether the system is metallic or insulating. However, 
since the conductance is given by a two-particle Greens 
function and involves complicated zero-temperature and 
zero-frequency limits it is difficult to calculate numeri- 
cally Q. 

In the following we concentrate on single-particle local- 
ization which is the most direct generalization of Ander- 
son localization to many-body systems. Experimentally, 
single-particle localization should be reflected in the tun- 
neling response of the system rather than in transport 
coefficients. 

The simplest measure of Anderson localization for a 
single-particle state |n) is the participation number P, 
defined as the inverse second moment of the spatial prob- 
ability distribution 

Pn^ElHj)! 4 (3) 
3 

where the sum runs over all sites j. In practice it is often 
averaged over all states with a certain energy e 

P-^) = ^^E^S(s-e n ). (4) 

A consistent generalization of this quantity to interacting 
systems should fulfill at least the following conditions: (i) 
it should be well defined for any many-particle state and 
(ii) it should reduce to (||) for non-interacting electrons. 
Moreover, the desired quantity should (iii) capture the 
physical idea of spatial localization and (iv) it should be 
easy to calculate. 

It has been suggested to define localization in a many- 
particle system via the spatial distribution of the charge 
difference between the states of the same many-particle 
system with N and N + 1 particles, respectively. While 
this quantity fulfills the above conditions (i), (ii), and 
(iv) it turns out that it is not a useful measure of local- 
ization in a disordered interacting system. The reason is 
that adding an extra electron to a disordered interact- 
ing system will very often not only add some charge at a 
few sites but completely rearrange the distribution of all 
electrons due to the frustration introduced by the com- 
petition between disorder and interaction. Thus simply 
calculating the participation number of the extra charge 
leads to an overestimation of derealization. This can 
already be seen at the example of the classical Coulomb 
glass where we know that the electrons are completely lo- 
calized. Nevertheless the ground states of systems with 
N and N + 1 particles can be drastically different so 
that the charge difference between the two is distributed 
on more than one site, effectively giving a participation 
number larger than one. In general, the method will al- 
ways fail, if adding an extra electron leads to a decrease 
of the charge at some particular site since then the charge 
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6 particles 



7 particles 



FIG. 1. Comparison of the ground states of quantum 
Coulomb glasses with 6 and 7 particles. The size of the circles 
is proportional to the charge density at the site. Note that 
adding an extra electron reduces the charge density at the site 
marked by an arrow. 



density difference is not a proper probability distribution 
anymore (see Fig. Q). 

There is, however, another quantity which fulfills con- 
ditions (i), (ii), and (iii) for a generalization of the inverse 
participation number. In particular, it nicely captures 
the physical idea of localization. This quantity is the 
probability R p for an electron to return to its starting 
site in infinite time. The energy-dependent return prob- 
ability can be expressed in terms of single-particle Greens 
functions 



p N ^s^o 7r n 



+ i5)G jj {e-i5). (5) 
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FIG. 2. Return probability for a quasiparticle in the quan- 
tum Coulomb glass and the Anderson model on a 3x4 lattice, 
t = 0.3, Wo = l,K = 0.5 



For non-interacting electrons P _1 (e) = R p (e). We note, 
however, that calculating R p requires the knowledge of all 
eigenstates of the Hamiltonian which makes this quantity 
numerically expensive. In Fig. ^ we show the numerically 
determined return probabilities of the quantum Coulomb 
glass and the Anderson model. The two systems behave 
very differently. For the Anderson model we obtain the 
well-known behavior of viz. a minimum in the 

band center and higher values in the tails. In the quan- 
tum Coulomb glass the return probability has a max- 
imum at the Fermi energy and decreases quickly with 
increasing distance from the Fermi energy. Thus we con- 
clude that in the quantum Coulomb glass single-particle 
excitations away from the Fermi energy tend to delocal- 
ize while the excitations close to the Fermi energy which 
dominate the low-temperature physics tend to localize 
We have investigated the values of R p close to the 
Fermi level for different values of the overlap t. In Fig. 
^ we show the resulting dependence and the correspond- 
ing data for the Anderson model. The figure shows that 
for t > 0.1 the interactions lead to stronger localization 
in agreement with the results of our Hartree-Fock cal- 
culation ||. For t = 0.1, quantum Coulomb glass and 
Anderson model show identical behavior within the nu- 
merical accuracy. Therefore, we cannot exclude that for 



c 
DC 



0.5 
0.4 

: 0.3 

L 

0.2 



□ 
♦ 



□ 


QCG 


♦ 


Anderson 



0.0 0.1 0.2 0.3 0.4 0.5 
t 

FIG. 3. Return probability for a quasiparticle at the Fermi 
energy as a function of t. Parameters are as in Fig. [| 
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FIG. 4. Fock space participation numbers for the ground 
states of the quantum Coulomb glass and the Anderson 
model. Parameters are as in Fig. 0. 



overlaps, and the authors found a delocalizing influence 
of the interactions. We note that, in principle, the con- 
cept of localization in Fock space can be applied to any 
problem in which the Hamiltonian can be decomposed 
into a reference part and a perturbation. Delocalization 
in Fock space then describes how the perturbation mixes 
the original eigenstates. In a recent letter [|l6| this con- 
cept has been applied to explain the transition in the 
width of excited states measured in tunneling conduc- 
tance experiments in quantum dots 

Let us further mention that the localization transitions 
discussed above also lead to transitions of the statistics of 
the corresponding energy levels [p|, p^|Jl^ , [l9[| . Analyzing 
the transitions of the level statistics is often numerically 
easier than dealing with the electron states itself. 



IV. CONCLUSIONS 



very small t, the behavior may turn around, i.e., the in- 
teractions may favor delocalization, see Ref. Jll| and the 
discussion of Fock space localization below. 

The return probability discussed above deals with the 
localization of a single quasiparticle. As already men- 
tioned, in an interacting system many-particle properties 
cannot in general be derived from single-particle ones. 
Therefore, to get a complete understanding of the local- 
ization properties one has to investigate many-particle 
quantities, too |15|. One approach which is in some 
sense complementary to the study of the return prob- 
ability is to analyze the localization properties of the 
many-particle states with respect to a Fock-space basis 
set {|a)}. Here the central quantity is the Fock-space 
participation number Pp of a many-body state \v) which 
is given by 



(6) 



If we chose the Fock space basis to consist of Slater de- 
terminants of site basis functions, |a) = cj . . . c\j|0) then 
Pp = 1 for completely localized electrons and Pp > 1 if 
the electrons can move. A measure like this has already 
been used to characterize the influence of long-range in- 
teractions on Anderson localization JjUfl . In Fig. [| we 
present our results for the Fock-space participation num- 
bers of the quantum Coulomb glass and the Anderson 
model. The data show the same tendency as the return 
probabilities at the Fermi level (Fig.|J) : For larger overlap 
t the interactions lead to stronger localization while for 
small t we cannot find a statistically significant difference 
between the quantum Coulomb glass and the Anderson 
model. Thus, a change in the behavior at very small 
t cannot be excluded. This would resolve the seeming 
disagreement with Ref. [jL5), where Fock space localiza- 
tion in a related model was investigated for very small 



In this paper we have discussed the generalization of 
Anderson localization to interacting systems. We first 
have considered the localization of single-particle excita- 
tions which can be described in terms of the return prob- 
ability R p . Experimentally, this quantity should be re- 
flected in the tunneling response of the system. We have 
also discussed the concept of Fock space localization for 
disordered interacting electron systems. As an example 
we have presented some numerical results for the quan- 
tum Coulomb glass model of disordered insulators. In 
this concluding section we will discuss some aspects that 
have not yet been covered. First, the quantum Coulomb 
glass is a model of spinless particles. In one-dimensional 
systems there is, however, some evidence that the elec- 
tron spin plays an important role in determining the be- 
havior of interacting disordered electrons ^o|. Therefore, 
in order to describe real electrons, including the spin into 
the quantum Coulomb glass will be necessary in the fu- 
ture. Second, our numerical examples were for very small 
lattice sizes. In order to quantitatively analyze the be- 
havior of disordered interacting electrons we have to ex- 
tend the calculations to larger lattices. This can either 
be done by making approximations such as Hartree-Fock 
Q or by developing better numerical algorithms. Third, 
one of the main tasks of the future will be to establish re- 
lations, if any, between the different types of localization 
and between the different quantities discussed here. 

This work was supported in part by the DFG under 
grant nos. Schr231/13-1, V6659/1-1 and SFB393 and by 
the NSF under grant no. DMR-95-10185. 



4 



[1] For reviews see e.g. P. A. Lee and T. V. Ramakrishnan, 

Rev. Mod. Phys. 57, 287 (1985), D. Belitz and T. R. 

Kirkpatrick, Rev. Mod. Phys. 66, 261 (1994). 
[2] For a recent review see B. Kramer and A. MacKinnon, 

Rep. Progr. Phys. 56, 1469 (1993). 
[3] M. Pollak, Discuss. Faraday Soc. 50, 13 (1970). 
[4] A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1975). 
[5] A. L. Efros, J. Phys. C 9, 2021 (1976). 
[6] Some recent results can be found in: A. Mobius, M. 

Richter and B. Drittler, Phys. Rev. B 45, 11568 (1992); 

T. Vojta, W. John and M. Schreiber, J. Phys. Condens. 

Matter 5, 4989 (1993); M. Sarvestani, M. Schreiber and 

T. Vojta, Rev. B 52, R3820 (1995). 
[7] M. Pollak, Phil. Mag. B65, 657 (1992). 
[8] See, e.g., K. M. Abkemeier, C. J. Adkins, R. Asal, and 

E. A. Davis, J. Phys. Condens. Matter 4, 9113 (1992); 

A. Aharony, Y. Zhang and M. P. Sarachik, Phys. Rev. 

Lett. 68, 3900 (1992); J. G. Massey and M. Lee, Phys. 

Rev. Lett. 75, 4266 (1995). 
[9] F. Epperlein, M Schreiber and T. Vojta, Phys. Rev. B 



(1997), in press, |cond-mat/9704068 
this volume. 



and contribution in 



[10] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[11] There is a relation between the conductance and the sen- 
sitivity of the ground state energy to twisted boundary 
conditions, see [Q, which can be used to calculate the 
conductance numerically. In dimensions larger than one, 
however, this method has its own problems. 

[12] D. J. Scalapino, S. R. White, S. C. Zhang, Phys. Rev. B 
47, 7995 (1993). 

[13] J. Talamantes, M. Pollak, and L. Elam, Europhys. Lett. 
35, 511 (1996). 

[14] Note that the opposite was suggested by I. L. Aleiner and 
B. I. Shklovskii, Int. J. Mod. Phys. B 8, 801 (1994). 

[15] There has been a very active line of research considering 
the localization of just two interacting electrons. While 
this is not a many-particle problem, it tries to capture 
essential aspects of the interplay of localization and in- 
teractions; see: D. L. Shepelyansky, Phys. Rev. Lett. 73, 
2607 (1994); Y. Imry, Europhys. Lett. 30, 405 (1995); 
D. Weinmann, A. Miiller-Groeling, J. L. Pichard, and K. 
Frahm, Phys. Rev. Lett. 75, 1598 (1995); F. v. Oppen, T. 
Wettig, and J. Miiller, Phys. Rev. Lett. 76, 491 (1996); 
R. A. Romer and M. Schreiber, Phys. Rev. Lett. 78, 515 
(1997). 

[16] B. L. Altshuler, Y. Gefen, A. Kamanev, and L. S. Levi- 
tov, Phys. Rev. Lett 78, 2803 (1997). 

[17] U. Sivan, F. P. Milliken, K. Milkove, S. Rishton, Y. Lee, 
J.M. Hong, V. Boegli, D. Kern, and M. DeFranza, Euro- 
phys. Lett. 25, 605 (1994). 

[18] B. I. Shklovskii, B. Shapiro, B. R. Sears, P. Lambrianides, 
and H. B. Shore, Phys. Rev. 47, 11487 (1993). 

[19] E. Hofstetter and M. Schreiber, Phys. Rev. B 48, 16979 
(1993). 

[20] T. Giamarchi and H. J. Schulz, Phys. Rev. B 37, 325 
(1988). 



5 



